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Abstract 

We study leading-order many-body effects of longitudinal optical (LO) 
phonons on electronic properties of one-dimensional quantum wire systems. 
We calculate the quasiparticle properties of a weakly polar one dimensional 
electron gas in the presence of both electron-phonon and electron-electron 
interactions. The leading-order dynamical screening approximation (GW 
approximation) is used to obtain the electron self-energy, the quasiparticle 
spectral function, and the quasiparticle damping rate in our calculation by 
treating electrons and phonons on an equal footing. Our theory includes 
effects (within the random phase approximation) of Fermi statistics, Lan- 
dau damping, plasmon-phonon mode coupling, phonon renormalization, dy- 
namical screening, and impurity scattering. In general, electron-electron and 
electron-phonon many-body renormalization effects are found to be nonmulti- 
plicative and nonadditive in our theoretical results for quasiparticle properties. 
PACS numbers:71.38.+i, 73.20.Mf, 71.10.+X 
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I. INTRODUCTION 



There has been a great deal of recent interest in ultranarrow confined semiconductor 
systems, called quantum wire structures, where the motion of the electrons is essentially 
restricted to be one dimensional (ID) because of their potential, both, in basic physics 
and for applications to novel device concepts PJ. Thus, both from the fundamental and 
applied physics viewpoints, there is interest in understanding the electronic properties of 
quasi-one-dimensional quantum wires. The quantum wire structures are usually fabricated 
by lithographic techniques 0, epitaxial growth on tilted superlattices ||, overgrowth of 
preetched V-groove patterns f§], stress patterning ||, and cleaved edge overgrowth method 
H which has been used to produce narrow quantum wires with atomic scale definition. 

In a doped polar semiconductor (in which most of the quantum wire structures are 
fabricated) free carriers couple to the longitudinal optical (LO)-phonons of the underlying 
lattice through the long-range polar Frohlich coupling. The carrier-LO phonon interaction 
leads to polaronic many-body renormalization of the single particle free carrier properties, 
e.g., polaronic effective mass renormalization, lowering of the effective band edge, inelastic 
scattering, broadening of the quasiparticle spectral function, etc. Even though the weakly 
polar III-V materials have rather small Frohlich coupling constants, the electronic properties 
can still be substantially modified by the Frohlich interaction, and, conversely, the phonon 
properties are affected by the Coulomb interaction between the electrons. Coulomb and 
Frohlich interactions cannot be disentangled in experimental situations, because many-body 
renormalization effects due to electron-electron and electron-phonon interactions in such a 
coupled system are fundamentally nonmultiplicative in nature (as well as being nonadditive). 

Our goal is to study a coupled ID electron-phonon many-body system treating electrons 
and phonons on an equal footing. Direct interaction between electrons via the Coulomb 
interaction and polar electron LO-phonon interaction (which is fundamentally Coulombic 
in origin, arising from the dynamical interaction between charge carriers and lattice ions) 
via the Frohlich coupling are among the most extensively studied many-body interactions 
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in solid-state physics. In a two-dimensional electron gas, quasiparticle properties of a polar 
electron gas have been extensively studied by treating the electron-electron and electron- 
phonon interactions on an equal footing |7J . In ref . |7| , Jalabert and Das Sarma find that the 
Coulomb and Frohlich interaction effects are nonmultiplicative, and the actual many-body 
correction for a 2D electron gas is substantially different from the one-polaron result (i.e. 
a single electron coupled to an LO-phonon system). In spite of substantial current interest 
in the properties of a polar quasi-one dimensional electron gas (ID EG) existing in GaAs 
quantum wires, there has been no detailed quantitative study of ID quasiparticle properties 
including both the electron-electron and electron-phonon interactions. In a GaAs (or, other 
III-V materials) quantum wire, the electronic energy scales (Fermi and plasmon energies) 
can be comparable to the LO-phonon energy, which produces strong mode coupling between 
plasmons and LO-phonon modes ||. Thus, a complete many-body analysis based on treating 
electrons, phonons, and plasmons equivalently within the same approximation scheme is 
needed. We provide such an analysis in this article based on the leading-order many-body 
perturbation theory. The Frohlich electron-LO phonon coupling constant in GaAs is small, 
making a weak coupling diagrammatic expansion for calculating quasiparticle properties a 
meaningful approximation. 

Well-established theoretical results PJiOj (Tomonaga-Luttinger model) predict unusual 
non-Fermi liquid-like properties for an interacting ID system. No matter how weak the 
electron-electron interaction the non-interacting Fermi surface is nonperturbatively unsta- 
ble in the presence of interaction. While the instability of a noninteracting ID Fermi gas to 
an interacting ID Tomonaga-Luttinger liquid is a well-established theoretical concept, the 
actual quantitative effect associated with the nonperturbative disappearance of the Fermi 
surface is, in fact, rather small at finite temperatures in a weakly interacting ID system such 
as a GaAs quantum wire ID EG. Most of the experimentally observed electronic properties 
of semiconductor quantum wires have so far been interpreted in the literature as effective ID 
Fermi liquids, albeit at finite temperature and in the presence of impurity scattering induced 
collision broadening. Indeed, based on a weak-coupling diagrammatic expansion in dynam- 
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ically screened Coulomb interaction, it has recently been argued |TT] that ID electrons in 



GaAs quantum wires behaves for most practical purposes as a thermally broadened Fermi 
liquid at finite temperatures and, most importantly, in the presence of impurity scattering 
effects which are invariably present in real semiconductor quantum wires. The physical 
reason behind this result is that impurity scattering cuts off the emission of arbitrarily low 
energy plasmons by overdamping them, and the instability of the Fermi surface, arising 
from the decay of quasiparticles into low energy plasmons, is prevented. Our approach in 
this paper is to include electron-LO phonon Frohlich interaction in this theory by doing an 
expansion in leading order dynamically screened Coulomb plus Frohlich interactions, both 
in the presence and absence of impurity scattering effects for the purpose of comparison. We 



note in this context that (1) the leading order calculation | llf , in the absence of impurity 
scattering, is consistent with the nonperturbative Tomonaga-Luttinger result in the sense 
that it indicates a non-Fermi liquid type behavior (with different analytic properties, how- 
ever); (2) the nonperturbative Tomonaga-Luttinger liquid type calculations have only been 
carried out for the point contact model electron-electron interaction whereas our work uses 
the true long range Coulomb (and Frohlich) interaction as the bare interaction. 

Usually the quantum wire confinement in the ^-direction (taken to be the growth direc- 
tion) is much stronger than the confinement in the y direction in real systems (assuming 
the 2D EG to be confined in the x — y plane). The electron is typically confined in the 
z-direction on the order of less than 100 A, whereas currently the best confinement in the y 
direction is approximately 300 A, leading to at least an order of magnitude difference in the 
energy-level spacings of the y (a few meV) and the z (the order of tens of meV) direction. In 
this paper we therefore assume for the sake of simplicity that the electron gas has zero thick- 
ness in the z direction, but has a finite width in the y direction. This, in fact, is one of the 
common models used in studying ID quantum wire structures. Confinement of the electrons 
in the y and z directions leads to the quantization of energy levels into different subbands. 
We assume the ID quantum limit where the energy separation between the lowest-energy 
(ground) and higher-energy (excited) subbands is so much larger than all other energy scales 



in the problem that the higher subbands can be ignored. In this paper we calculate in the ID 
quantum limit the zero temperature electron self-energy from the effective total dynamical 
electron-electron interaction, which is calculated by taking into account both the Coulomb 
electron-electron interaction and the LO phonon-mediated Frohlich interaction, within the 
random-phase approximation (RPA). We treat electron-electron and electron-phonon in- 
teractions equally, and our only approximations are the use of the RPA and the neglect of 
vertex corrections, i.e. we keep the lowest order diagram in the total effective dynamically 
screened interaction. This approximation has long been employed successfully to calculate 
properties of three- and two-dimensional electron systems fl2yi3f. Furthermore, the RPA 



gives the exact low-energy screening properties of one- dimensional systems [J4| . The Feyn- 
man diagrams for the self-energy and the effective interaction in these approximations are 
shown in Fig. 1. We ignore diagrams which are higher order in the screened interaction, 
such as the diagram shown in Fig. 1(c). 

This paper is organized as follows. In Sec. II we introduce our model for the coupled ID 
electron-LO phonon system. In Sec. Ill we describe our self-energy and spectral function 
calculations for the quantum wire at T = 0. In Sec. IV we present our numerical results for 
various quasiparticle properties. We provide a conclusion in Sec. V. Throughout this paper 
we have used the terminology in which we refer to the situation having both Coulomb and 
Frohlich interactions as the coupled system, whereas the situation with only the Coulomb 
interaction (i.e. with the Frohlich coupling turned off) is called the uncoupled system. Also 
note that in this paper, h = 1. 



II. MODEL 

Our model consists of a ID EG coupled to bulk dispersionless LO-phonons at zero temper- 
ature. Electrons interact among themselves through the Coulomb interaction and through 
virtual-LO-phonon exchange via the Frohlich interaction. In calculating the effective ID 
electron-phonon interaction we sum over the phonon wave vector in the other two dimen- 
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sions in the standard manner 0. The Coulomb interaction is logarithmically divergent in 
the ID wave vector space, and therefore we use the more realistic finite width quantum 
wire model ]nj. In our extreme quantum limit model where only the lowest ID subband 
is occupied by the electrons, we obtain the Coulomb interaction matrix element V c (q) by 
taking the quantizing confinement potential to be of infinite square well type [n],rf|, 



2e f 1 

V c (q) = / dx K (qax) 

Coo Jo 



[l - x) (2 + cos(27ra)) + — sin(27rx) 



(1) 



where K the zeroth-order modified Bessel function of the second kind, a the width of the 
ID EG (i.e. the effective confinement width), and is the high frequency background 
lattice dielectric constant. In the long wavelength limit its asymptotic form becomes 



V c (q) = —K (qa). 



(2) 



The LO-phonon mediated electron-electron interaction is dependent on wave vector and 
frequency, 



V ph (q,u;)=MlD (uj) 



M q is the effective ID Frohlich matrix element given by 



(3) 



r 2 , / \ ^LO 



Mt = V c (q)- 



1 



(4) 



where c^lo is the LO-phonon frequency, eo is the static lattice dielectric constant. The 
unperturbed retarded bare LO-phonon propagator is given by 

2u; LO 



(5) 



uj 2 — cg>lo + i 0+ sgn(u;) ' 
The effective electron-electron interaction is obtained in the RPA by summing (Fig. 1(b)) 
all the bare bubble diagrams, 

V c (q)+V ph (q,ou) 



V eS (q,uj) 



l-[V c {q)+V ph (q,cu)]U (q,u;) 
et(?,w)' 



(6) 



where Ho(q, to) is the complex irreducible ID polarizability given by the bare bubble diagram. 
The analytic form of Ho(q, to) for complex frequency is given by |T5| 

z 2 — (q 2 /2m — qvF) 2 



m 



U (q,z) = — In 
Tiq 



(7) 



z 2 — (q 2 /2m + qvp) 2 

where the principal value of logarithm [i.e. —n < Im[ln(x)] < 7r) should be taken. In evalu- 
ating U Q (q,uj) for real frequency, the limit z = uo + i0 + is taken, (vp is the Fermi velocity). 
€t(q,uj) is the total dielectric function, which has an electron and a phonon component, 

1 - eoo/eo 



e t (q,u) = 1 - V c (q)U (q,u) + 



(8) 



eoo/eo - u 2 / u)l - i sgn(w) 0+ ' 

We include collisional broadening effects arising from electron- impurity scattering by ap- 
propriately modifying |lTJ our dielectric function. Impurity effects in screening are usually 
introduced diagrammatically by including impurity ladder diagrams into the electron-hole 
bubble. Since the exact expression for the polarizability within this ladder vertex diagram 



approach is complicated, we use the particle-conserving polarizability function, H 



given by Mermin |T7|]. Impurity scattering self-energy effects are not included in the single 
electron Green's function. In the presence of an impurity scattering induced level broadening 
of 7, the polarizability function is given by, 

(u + ry)n (g,u; + ij) 



w + i 7 [n (g,u; + 27)/no(g, o)]' 



(9) 



III. ELECTRON SELF-ENERGY 

The electron Green's function G(k, u), or equivalently, the electron self-energy, E(A;, u) = 
Gq x (k, uj)—G~ l {k 1 u>), where Go(k, to) is the bare noninteracting Green's function, provides a 
complete description of all single electron properties of an interacting many-electron system. 
Knowledge of G(k,u) or E(/c, u) allows one to calculate many experimentally observable 
one-electron properties of a system, such as the spectral density function, effective mass, 
scattering rate, lifetime, distribution function, band gap renormalization, etc. The electron 
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self-energy of the coupled system within the leading-order (Fig. 1(a)) GW approximation 
fl8| neglecting vertex corrections at T = is given by 



E*(ife,u;) = i 



30 dq f°° dJ . t 

00 Z7T J-00 Z7T 



(10) 



where G l is the Green's function for the noninteracting electron gas (the superscripts t are 
to denote time-ordering; otherwise, the functions in this paper are retarded) 

1 



& {k,u>) 



(11) 



to — £(k) + z0+sgn(/c — hp) ' 

with = (A; 2 — hp) /2m. The time-ordered quantities are easily related to the retarded 
quantities pTOf] . The self-energy can be separated into a frequency-independent exchange 
contribution and a correlation part, 



E(&, lo) = E cx (/c) + E cor (A;, uj). 

The exchange or Hartree-Fock contribution is given by 

dq 



(12) 



E«(*0 = - / r^ F (A; + g)V c (g), 



(13) 



where np(/c + g) = 0(A;j? — |fc + g|) is the Fermi function at T = 0, and E cor (fc,u;) is defined 
to be the part of H(k,u) not included in E ex (A;). In the GW approximation, the E cor (A;, u;) 



can be written in the line and pole decomposition [10,19 



E cor (A;, u) = Si ine (/c, to) + E pole (A;, u) 



(14) 



where 



Ei ine (/c, u) 



30 dq^ . r°° duj' 1 

00 2ir ° J -00 2tc + g) — w) — iui' 



- 1 



(15) 



and 



E pole (A;, u) = r p- [6 (u - £(* + g)) - 5 (-£(* + g))] 



£t(q,£(k + q) -u) 



- 1 



(16) 
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The Sii ne (/c, a;) is completely real because e t (q, iuj') is real and even with respect to uj' . Thus, 
the imaginary part of the self-energy comes entirely from Im[S pole ]. The pole contribution 
of the self-energy can be written as 



Spole(^^) 



fco M 



e(u>) 



9{-uj - E F ) 



6{uj + E F )0(-u) 



F(q,k,uj)dq, (17) 



where k (uj) = ^2m(uj + E F ) and F(q,k,uj) = ^[f(q,k,ui) + f(—q,k,ui)] where f(q,k,uj) 
is given by 



f(q,k,uj) = V c (q-k) 



- 1 



(18) 



_€ t (q-k,£(q) -uj) 

Quasiparticle excitation energies uj q of the system are obtained from the real part of the 
self-energy Re[S(/c,cj)] by solving Dyson's equation, 



,(Jfe) -f(fc) - ReE(k,uj) = 0. 



(19) 



Once the self energy H(k,u) is known the single particle spectral function A(k,u>) is readily 
calculated. A(k, uj) contains a important information about the dynamical structure of the 
system and is given by 

2\lmE{k,uj)\ 



A{k,uj) 



{uj - £(k) - ReS(/c, uj)} 2 + {Im£(fc, uj)} 2 ' 



(20) 



If Im[S(/c,a;)] goes to zero faster than as uj — > uj^k), the spectral function becomes a 
(^-function peak at the quasiparticle energy uj q , 



A{k,uj) = Z(k) 5{uj -uj q (k)). 



(21) 



The strength of the peak is given by the quasi-particle renormalization factor, 

n -l 



Z{k) 



( _ dReJ:(k,uj) i 
du 



(22) 



The characteristic of a Fermi system is that it has a finite discontinuity at the Fermi momen- 
tum in the momentum distribution function. If Z(k F ) ^ 0, the electron distribution as a 
function of momentum has a finite discontinuity at the Fermi momentum. Thus, if A{k F , a;) 



shows a ^-function peak at k — k,F and uj = the system is a Fermi liquid. However, if 
Im[E(/cF,u;)] goes to zero slower than uj as uj — > 0, then the spectral function = 0) 

has finite broadening at and is not a 5-function, implying that the system is not a Fermi 
liquid. 

The spectral function A(k, uj) can be interpreted as a probability that an electron has a 
momentum k and energy uj, and satisfies the sum rule 



which is equivalent to the conservation of the number of particles (electrons) as the interac- 
tion is turned on, starting from a non-interacting picture. In all our calculations, this sum 
rule is explicitly checked and is found to be satisfied to within less than a percent. The 
quasiparticle broadening or the damping rate T(k) is given by the imaginary part of the 
self-energy, 



where = k 2 /2m — /i is the quasiparticle energy measured with respect to the chemical 
potential //. From the damping rate we can calculate the quasiparticle scattering rate 2T(k), 
the inelastic life time r(k) = [2T(k)]~ 1 , and the inelastic mean free path l(k) = v(k)r(k), 
where v(k) is the electron velocity. 



In Fig. 2 we show the real and the imaginary parts of the self-energy and spectral 
functions for uncoupled (thin lines) and coupled (thick lines) systems without any impurity 
scattering. We show the self energies and spectral functions only at k — (band edge) and 
k = kF (Fermi energy). We use Ep [kp) as our unit of energy (momentum) throughout our 
discussions. In this calculation we use the parameters corresponding to GaAs: m = 0.07m e 
(m e is the free electron mass), e = 12.9, = 10.9, and u LO = 36.8 meV. The well width of 




(23) 



T(k) 



■Im£(fc,£(fc)), 



(24) 



IV. NUMERICAL RESULTS 
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a = 100A and the electron gas density of n = 0.56 x 10 6 cm -1 , which corresponds to a Fermi 
energy Ep ~ 4.4meV and a dimensionless density parameter r s = 4me 2 /nkpeo = 1.4 with 
fcp = im/2, are used for this calculation. The intersections of Re[E(/c, uj)] and the straight 
line uj — indicate the solutions to Dyson's equation. For k = and the uncoupled 
system we find three solutions of the Dyson's equation (Eq. 18). Two solutions (the first 
and the third one from u = 0) show the 5-function peaks in the spectral function, but the 
second solution has a broad incoherent structure indicating finite damping in the spectral 
function. The first 5-function peak in A(k,u) (with Z = 2tt x 0.33), which corresponds to 
the usual quasi-particle (i.e., a bare particle surrounded by a cloud of virtual plasmons and 
particle- hole excitations), is slightly shifted from the noninteracting single particle energy 
uj = The second 5-function peak (with Z = 2tt x 0.31) has been called a plasmaron 

which is interpreted as a hole coupled to a cloud of real plasmons. For the Fermi wave vector 
(k = kp) there is a only one solution to Dyson's equation at uj = and a strong peak in 
A(k,uj). However, this peak is not the 5-function peak, since Im[£(fcp, uj)] — > uyhiu as 
uj —>■ | Thus the uncoupled ID EG without any impurity scattering is not a Fermi 



liquid even within the leading-order GW approximation [11], and does not have well-defined 



quasiparticles. This was already reported in ref. [11 



Since lm[e c (q, uj)}^ 1 ^ (the subscript c stands for the Coulomb part; i.e., the uncoupled 
system) within the electron-hole continuum, the contribution to the imaginary part of the 
self-energy, Im[E(fe, uj)] = Im[E po i e (A;, a;)], comes from both the electron-hole continuum and 
the collective plasmon mode. Im[e c (g, uj)]^ 1 can be written as 
1 -e 2 



Im 



e c (g,w) 



f 2 4- f 2 



[0(|u;i| - uj)6(uj - \uj 2 \)9(uj) + 6(-\uj 2 \ - uj)6{uj + \ui\)e(-u)] (25) 



—71 



du 



5{uj - u p i), 



where ei = Re[e(g, uj)] and e 2 = lm[e(q,uj)], and the plasmon dispersion, uj p \, is given by [15 



^pi(?) 



(26) 



g(q) - 1 

where g(q) = exp[q7i/mV c ], U\ = 2q + q 2 , and uj 2 = 2q — q 2 . The first (second) term 
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in Eq. (25) gives the electron- hole (plasmon) contribution. From Eq. fll7|), we see that 
the non-zero electron-hole contribution shows up in Im[£(/c, uj)] only if uj > 1 — k 2 for 
uj > and if — (k + 1)(A; + 3) < uj < — 11 — k 2 \ for uj < 0. The plasmon part (the 
second term in Eq. (25) contributes to Im[E(A;, a;)] if uj > uj p \(1 ± k) for uj > and if 
— uj p \(1 ± k) < uj < — 11 — k 2 \ for a; < 0. Thus, for k = 0, Im[E(fc, u;)] 7^ in the range 
— uj p \(1) < uj < — 1 and u; > 1 (the electron-hole contribution, uj > 1 and — 3 < uj < — 1, 
and the plasmon contribution, — cj p i(1) < uj < — 1 and a; > co> p i(l)). Since the plasmon 
contributions turn on at uj = ±u; p i(l) the finite discontinuities in Im[E(fc, uj)] take place at 
those points with the magnitude ir/^\u=±u pl . A finite discontinuity in Im[S] gives rise to a 
logarithmic singularity in Re[S], which can be verified using the Kramers- Kronig relation. 
For k = kp-, Imp] 7^ in the range uj > and —uj p i(2) < uj < ( the electron-hole 
contribution, uj > and — 6 < uj < 0, and the plasmon contribution, uj > u> p \(2) and 
— cj p i(2) < uj < 0). Turning on the plasmon contribution at uj = ±cj p i(2) induces the finite 
discontinuity in Im[S] and a logarithmic singularity in Re[S]. 

For the coupled ID EG and k = we find three undamped solutions and two damped 
solutions in Dyson's equation. The first 5-function peak in A(k, uj) near uj = 0, the quasi- 
particle peak, occurs at almost the same energy as the uncoupled case, but with slightly 
higher strength (Z = 2ir x 0.36). The second 5-function peak, the plasmaron peak, is 
shifted toward the Fermi energy and has much weaker strength (Z = 2n x 0.22) than in the 
uncoupled ID EG, since the coupled plasmon-like mode has a smaller energy and weaker 
strength than in the uncoupled case ||. Unlike the uncoupled system we now find a third 
5-function peak with strength Z = 2ir x 0.05 near the frequency uj = —Ep — c^lo This 
phonon satellite can be interpreted as being a "phononron" similar to the plasmaron, that 
is, a quasi- hole coupled to real LO-phonons. For k = kp we find only two peaks : the 
quasi-particle peak at uj = and the phonon satellite (the 5-function peak with strength 
Z = 2n x 0.03). As uj —>■ the main contribution to Im[S] comes from a plasmon, since 
the phonon part has a gap at q = and the integration is performed near q = 0. Thus, the 
behavior of Im[E] for the coupled system is the same as that for the uncoupled system and 
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the quasi-particle like peak near oj = is not a strict 5-function peak, which means that 
the coupled ID EG with electron-phonon interaction is not a Fermi liquid either, and has 
no true long-lived quasi-particles. 

For the coupled system there is also a direct phonon contribution to the imaginary part 
of the self-energy along with the electron-phonon and plasmon contribution discussed above. 
Im[e t ] _1 can be written as 
1 -e 2 



Im 



where and are the plasmon-like and the phonon-like collective modes in the coupled 
system, respectively |§. Their behaviors in the long wave length limit (q — > 0) are given by 



4 + 4 



[0(|wi| - oj)6(uj - \u 2 \)9(u) + 6(-\u 2 \ - oj)6(oj + |u;i|)0(-u;)] (27) 



— IT 



,(+) 



duj 



,(-) 



-i 



8(u - U^) - IT 



duj 



,(+) 



-i 



i + 



OJ 



LO 



TO 



2r s \\n(qa)\q 2 + 0(q 4 



LO 



(28) 



and 



UTO t 



2qJr s \\n(qa)\+0(q 3 



(29) 



The contribution to Im[£] from the electron-hole continuum is the same as that for the 
uncoupled system, but the contribution of the plasmon-like mode, C0p[\ is the same only 
as long as to < u)^ 0) since the plasmon-like mode ceases for to > c^to [See Ref. [|J]. The 
discontinuities in Im[E] induced by this mode move toward the Fermi energy (uj = 0) and the 
strength of the step becomes smaller, since |o; p i| > l^ij |. The phonon-like mode contributes 



to Im[S] if the following conditions are satisfied: if u > 0, q 2 



LU + OJ. 



(+) 
pl 



with 



1 < q < \/u + 1, and if uo < 0, q 2 — 1 — oj — oj^' = with < q < 1. From these 
conditions we find the range in which the phonon-like mode contributes to the Im[S] to be 
oj > cjpi (1 ± k) and —1 — oj^ \k) < uo < —uo^\l — k). As for the plasmon, turning on the 
phonon-like mode produces a discontinuity in Im[E] and a singularity in Re[E]. For k = kp, 
since ^^(1 ± k) > u>ro, the abrupt steps in Im[£] for oj > are due to the phonon-like 
mode. 
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Fig. 3 shows the results for the coupled ID EG with impurity scattering effects (7 = 
Ep). The inclusion of impurity scattering introduces collisional damping in the collective 
modes. Since the ID plasmon is strongly affected by impurity scattering, (i.e., they become 
overdamped at small q), the damping of the plasmaron (the second peak in A from uo = 0) 
peak is stronger than that for other excitations. However, the phonon part of the dynamical 
dielectric function is not affected by the impurity scattering. Thus, the phonon satellite (the 
third peak from uo = in A) still shows a sharp peak ( almost undamped ), and the abrupt 
steps in Imp] still exist at uo = ±co>lo and uo = —Ep — culo- The discontinuities in Re[E] 
and the singularities in Re[S] induced by the plasmon-like mode are both suppressed by the 
impurity scattering. For k = kp, the contribution to Im[£(A; = kp,uo)] at u = mainly 
comes from the single electron-hole pair excitations as in the uncoupled ID EG, since the 
long wavelength plasmon is now overdamped and the phonon has a gap at a small q. As 
uo — > the contribution to Im[E] from the single-particle excitations is given by uo 2 \ lnu;| 3 , 
and A(kp,u) ~ |ln(|u;|)| 3 + 2tiZ f 5{uo) [|TTJ. Our direct numerical calculation shows that 
the weights of the S function are 2ir x 0.31 for the uncoupled system and 2tt x 0.33 for the 
coupled system. Therefore, the coupled ID EG with impurity scattering is a Fermi liquid 



similar to the behavior found in ref. [[□]] for the uncoupled system in the presence of impurity 
scattering. 

In Fig. 4 we show the strength Z(k) of the undamped excitations for uncoupled (thin 
lines) and coupled (thick lines) system without impurity scattering as a function of wave 
vector. As the wave vector increases the strength of the ^-function peaks for the plasmarons 
(dashed lines) decrease quickly and vanish at k — 0.75k F for the uncoupled system and 
at k = 0.3kp for the coupled system, but the strengths of the regular quasiparticle peaks 
(solid lines) and the phononron peak (dot-dashed line) decrease slowly. On the other hand, 
as the wave vector approaches the Fermi momentum the strength of the quasiparticle peak 
decreases and vanishes at the Fermi momentum. Analytically we can show that this behavior 
is logarithmic as in the uncoupled system. The most dominant term in the derivative of the 
real part of the self-energy comes from E po i e as u — > 0. From Eq. (16) we have as uo — ► 
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dq5(u - £ k+q )v c (g) 



1 

1 



«c(S 1-7^ - l) + «o(2) - 1 



2 y Ut(f,0) J cv/ ie 4 (2,0) 
Thus, as u — > the first term shows a logarithmic singularity. 



(30) 



k=kF ~ - (In M + 1) . (31) 



duj 

This demonstrates that even within our leading order approximation the coupled system 
without impurity scattering is not a Fermi liquid, since the renormalization constant at 
k = kp vanishes, (i.e., the momentum distribution is continuous at the Fermi momentum). 

Fig. 5(a) shows the quasiparticle damping rate Y{k) for both the coupled system 
(thick lines) and the uncoupled system (thin lines) with the parameters a = 100A and 
n = 0.56 x 10 6 cm _1 . (Inset in Fig. 5(a) shows the plasmon-phonon mode coupling.) The 
corresponding inelastic mean free paths Ik are shown in Fig 5(b). For the uncoupled ID EG 
the quasiparticle scatters by plasmon emission, which corresponds to the sharp threshold 
in the thin solid line shown in Fig. 5(a). There is no ID single particle electron-electron 
scattering below a threshold critical wave vector because the conservation of energy and 
momentum restricts electron-electron scattering to an exchange of particles. In 2D and 
3D systems there is allowed scattering below the threshold critical wave vector due to the 
excitation of electron-hole pairs. For the coupled ID EG the quasiparticle decays via the 
emission of coupled plasmon-phonon mode |20| (u;_ and uo + in the inset of Fig. 5(a)), which 
corresponds to the two peaks in the thick solid line of Fig. 5(a). The first step corresponds 
to the a;_-emission threshold and is located below the k c of the uncoupled ID EG because 
uj- < ujq. The second step corresponds to the co> + -emission threshold and occurs at a wave 



vector larger than k = \/ujlo + 1 because lu + > ujlo- When only Frohlich interaction in the 



ID quantum wire is considered [21 1 the LO-phonon emission threshold has been predicted 



at the fixed critical wave vector k = \/ujlo + 1- But in our coupled model the locations 
of C(j-i--emission threshold depend on the density of the electron as in the plasmon emission 
threshold. The thin lines in Fig. 5(a) show the damping rate with the impurity scattering 

15 



effects (7 = Ep). Since the plasmon line is broadened by impurity scattering there are 
inelastic scattering events below the critical wave vectors. The breaking of translational in- 
variance by impurities, furthermore, relaxes momentum conservation and permits inelastic 
scattering via single-particle excitation (except at the Fermi surface where the quasiparticle 
is always undamped). Thus, the nature of the sharp threshold in our doped system is totally 



different from that discussed in ref. which considers undoped quantum wires. The sharp 



divergent thresholds obtained in Ref. |21| result from the divergence in the density of states 
at the bottom of the band of one- dimensional systems, and therefore in doped samples, when 
the bottom of the band is no longer accessible, these sharp thresholds do not occur. In our 
theory, the divergences are due to the divergence in the joint density of states at the coupled 
plasmon-phonon emission threshold, and they occur at any doping. 



V. SUMMARY AND CONCLUSION 

In this paper we have calculated the quasiparticle self-energy, the spectral function, and 
the damping rate within the leading order dynamical screening approximation by treat- 
ing electron-electron and electron-phonon interactions on an equal footing both with and 
without impurity scattering effects. Our many-body theory includes the important physical 
effects of dynamical screening, phonon self-energy correction, plasmon-phonon mode cou- 
pling, Fermi statistics, and Landau damping. Even though the problem is treated within the 
GW framework of the leading order effective interaction approximation (i.e. RPA screening 
and neglect of vertex corrections), our results should be quite valid in GaAs-based ID EG. 
This is because GaAs has a very weak Frohlich coupling, which justifies the neglect of the 
electron-phonon vertex corrections, and its low effective mass and the large dielectric con- 
stants gives a small effective r s parameter, (and in this limit generally RPA is valid), which 
makes direct electron-electron interaction also weak in the perturbative sense. 

The quasiparticle properties at uj ~ are affected only by low energy processes and 
hence are not changed much by the inclusion of the optic-phonon coupling because the 
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optic-phonon modes have a gap to excitation. The optic-phonon coupling does, however, 
shift some of the spectral weight to an extra undamped excitation at uj ~ —(E F + u LO ), 
which is the phonon analogy ( "phononron" ) to the plasmaron excitation, and it shifts the 
plasmaron energy up slightly, due to level repulsion between the plasmaron and the phonon 
excitation. As in the uncoupled system, inclusion of impurity scattering tends to smear out 
the above features. 

The properties at the Fermi surface, such as the impurity-driven restoration of the dis- 
continuity in the distribution function, are virtually the same for the coupled and uncoupled 
systems, again because of the gap of the optic-phonon mode. At higher energies, the coupled 
system has an additional phonon-like mode oj + which can scatter electrons. This leads to 
another divergence in the scattering rate Y{k) (i.e., in addition to the one from plasmon 
emission). Note that this purely many-body divergence occurs only when one couples the 
electron-electron and electron-phonon interactions together; no such divergence is obtained 
for calculations in doped systems using only the Frohlich interaction. 
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FIGURES 

FIG. 1. (a) Electron self-energy in leading order in the effective dynamical interaction, (b) 
Effective dynamical interaction (thick wiggly line) V e s calculated in the RPA. Thin wiggly lines 
(dashed lines) represent the Coulomb electron-electron interaction V c (the LO-phonon-mediated 
electron-electron interaction Vph), and the bubble the irreducible polalizability IIo. (c) Higher-order 
self-energy diagrams neglected in our calculation. 

FIG. 2. Self energy S(fe, uj) ( (a) for k = and (b) for k = kp) and spectral function A(k, uj) ( 

(c) for k = and (d) for k = &f ) as functions of the frequency uj without any impurity scattering 
(7 = 0). The straight lines are given by uj — £k — A*, and their interactions with Re[S] indicate the 
solutions to Dyson's equation and correspond to a quasiparticle peak. For clarity we plot |ImS|. 
Thin (thick) lines correspond to the uncoupled (coupled) ID EG. The vertical arrows in (c) and 

(d) denote 5 functions in the spectral function. For k = 0, we find at uj = -4.85 and -0.96 with 
weights 2ir x 0.31 and 2ir x 0.33 respectively for the uncoupled system, and at u = -9.65, -4.28, 
and -0.96 with weights 2ir x 0.05, 2tt x 0.22, and 2tt x 0.36 respectively for the coupled system. 
For k = kF, we find a peak at w = —9.81 with weight 2ir x 0.03 only for the coupled system. The 
parameters corresponding to GaAs are used: m = 0.07m e (m e is the free electron mass), eo = 12.9, 
£oo = 10.9, and wlo = 36.8 meV. The well width of a = 100A and the electron gas density of 
n = 0.56 x 10 6 cm _1 , which corresponds to a Fermi energy Ep m 4.4meV and a dimensionless 
density parameter r s = 4me 2 /irkpeo = 1.4 with fc^ = irn/2, are used for this calculation. 

FIG. 3. Self energy T,(k,oj) ( (a) for k = and (b) for k = kp ) and spectral function A(k,u>) 
( (c) for k = and (d) for k = kp ) as functions of the frequency uj with impurity effects (7 = Ep). 
Thin (thick) lines correspond to the uncoupled (coupled) ID EG. For k = kp, vertical arrows in 
the spectral function at uj = denote 5 functions with weight 2ir x 0.31 for uncoupled system and 
2ir x 0.33 for uncoupled system. 
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FIG. 4. Strength of the undamped excitations (the quasiparticle renormalization factor Z(k)) 
for uncoupled (thin lines) and coupled (thick lines) system without impurity scattering as a function 
of wave vector. Solid lines (dashed lines) denote quasiparticle peaks (plasmarons) and dot-dashed 
line a phononron. Inset shows the logarithmic approach of the renormalization factor to zero as 
k — > kp. 

FIG. 5. (a) Damping rate T(k) and (b) the corresponding mean free path l(k) = v(k)/2T(k) 
as a function of k. Solid (dashed) lines indicate the coupled (uncoupled) system and thick (thin) 
lines for 7 = (7 = Ep). The inset in (a) shows the plasmon-phonon mode coupling. 
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